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^yQ I Abstract: After inserting the heavy quark mass dependence into QCD partonic evolution equations, 

•^ ■ we determine the mean charged hadron multiplicity and second multiplicity correlators of jets produced 

in high energy collisions. We thereby extend the so-called dead cone effect to the phenomenology of 
Q ■ multiparticle production in QCD jets and find that the average multiplicity of heavy-quark initiated jets 

O . decreases significantly as compared to the massless case, even taking into account the weak decay prod- 

ucts of the leading primary quark. We emphasize the relevance of our study as a complementary check 
of 6-tagging techniques at hadron colliders like the Tevatron and the LHC. 
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1 Introduction 

Since the very beginning of cosmic and accelerator physics, the study of jets has been playing a promi- 
nent role in the rise and development of the Standard Model (SM) lU. For example, the observation 
of three-jet events in electron-positron collisions at DESY provided a direct experimental evidence of 
the existence of gluons. Nowadays, QCD furnishes the theoretical framework for jet analysis, and 
conversely, jet studies furnish precise tests of both perturbative and non-perturbative QCD, as well as 
constraints and determinations of QCD parameters. 

While high-energy hadronic interactions are dominated by the production of secondaries with rather low 
transverse momentum (j)t) with respect to the beam axis, high-pt jets are expected to become one of 
the cleanest signatures for New Physics (NP) to be discovered at the LHC. On the other hand, QCD 
processes often are the most important background of such NP signatures, therefore requiring a good 
understanding of QCD jet rates and features. 

High-pi jets can be initiated either in a short-distance interaction among partons of the colliding hadrons, 
or via electroweak (or new physics) processes. One well-known example is given by the decay chain of 
the top quark t — )• H^ b, where the b quark should start a jet. Thus the ability to identify jets from 
the fragmentation and hadronization of b quarks becomes very important for such Higgs boson searches. 
Needless to say, the relevance of 6-tagging extends over many other channels in the quest for new physics 
at hadron colliders. The experimental identification of 6-jets relies upon several of their properties in 
order to reject background, e.g. jets initiated by lighter quarks or gluons. First, the fragmentation is hard 
and the leading 6-hadron retains a large part of the original b quark momentum. In addition, the weak 
decay products may have a large transverse momentum with respect to the jet axis therefore allowing 
separation from the rest of the cascade particles. Lastly, the relatively long lifetime of ^-hadrons leading 
to displaced vertices which can be identified by using well-known impact parameter techniques tlj. Still, 
a fraction of light jets could be mis-identified as 6-jets, especially at large transverse momentum of the 
jet. Now, let us point out that an essential difference between heavy and light quark jets results from 
kinematics constraints: the gluon radiation off a quark of mass m and energy £^ ^> m is suppressed 
inside a forward cone with an opening angle 0^ = m/E, the so-called dead cone phenomenon EIH. 

In this paper, we compute the average (charged) multiplicity and multiplicity fluctuations of a jet initiated 
by a heavy quark. For this purpose, we extend the modified leading logarithmic approximation (MLLA) 
evolution equations ||5l to the case where the jet is initiated by a heavy (charm, bottom) quark. The 
average multiplicity of light quarks jets produced in high energy collisions can be written as N{Y) oc 
exp < / j{y)dy > (Y and y are defined in section|2ll, where 7 ^ 70 + A7 is the anomalous dimension 
that accounts for soft and coUinear gluons in the double logarithmic approximation (DLA) 70 ~ y/a^, in 
addition to hard coUinear gluons A7 ~ Os or single logarithms (SLs), which better account for energy 
conservation and the running of the coupling constant as [5]. Thus, in the MLLA both contributions can 
be added in the form 

7~ V^ + a, + 0(a3/2), ^y^ 

such that SLs corrections with respect to DLA are of relative order 0{y/a2). Notice that in jet calculus 
the scheme of resummation differs from that used in NLO hard scattering cross-sections calculations 
from the factorization theorem fSHHl, where no resummation over the soft gluon logarithms is performed. 



As a next step in our approach, we include the dead cone phenomenon into the massless quark equations 
by using the massive leading order (LO) splitting functions, first computed in QED for electrons and 
photons [9J, and by replacing the massless quark propagator l/k^ by the massive one 1/(A;^ + m?), as it 
was carried out for the evaluation of jet rates in the e+e~ annihilation ifTOl . Furthermore, we demonstrate 
that the whole phase space of the heavy quark dipole QQ produced in the e+e^ annihilation L4J reduces 
to that of the heavy quark jet event in the soft and coUinear limits. 

In the present work we evaluate the mean multiplicity and the second multiplicity correlator as a function 
of the mass of the heavy quark. Actually, the mass depends on the scale Q of the hard process in which 
the heavy quark participates. In the MS renormalization scheme, for example, the running mass becomes 
a function of the strong coupling constant as{Q), and the pole mass defined through the renormalized 
heavy quark propagator (for a review see ifTTl and references therein). Nonetheless, in our QCD leading- 
order analysis no distinction between running and pole mass of heavy quarks has to be made. On the 
other hand, the quark mass considered as a parameter in our calculations should not play a significant role 
in the jet evolution until the virtuality Q drops down to values of the order of the quark mass. Therefore, 
the dominance of the dead cone phenomenon at low Q justifies neglecting the running of the heavy quark 
mass, hence adopting a value close to the pole mass as a natural cut-off for final-state soft and collinear 
singularities. In order to assess the scheme-dependence of our calculations beyond leading-order we 
have considered a broad range of possible values for the charm and bottom masses: rric = 1 — 1.5 GeV 
and mf, = 3 — 5 GeV, respectively ifTTTl . 

Lastly, we will see that under the assumption of local parton hadron duality (LPHD) as hadronization 
model II121I13I . light- and heavy-quark initiated jets show significant differences regarding particle mul- 
tiplicities as a consequence of soft gluon suppression inside the dead cone. Such differences could be 
exploited by using auxiliary criteria complementing 5-tagging procedures to be applied to jets with very 
large transverse momentum, as advocated in this paper. 

2 Kinematics and variables 

As known from jet calculus for light quarks, the evolution time parameter determining the structure of 
the parton branching of the primary gluon is given by (for a review see |5 1 and references therein) 

y = \n(^\ k± = zQ>Qo, Q = Ee>Qo, (2) 

where k± is the transverse momentum of the gluon emitted off the light quark, Q is the virtuality of the 
jet (or jet hardness), E the energy of the leading parton, Qo/E < G < Go is the emission angle of the 
gluon (Q <C 1), Go the total half opening angle of the jet being fixed by experimental requirements, and 
Qq is the collinear cut-off parameter. Let us define in this context the variable Y as 

y = Y + lnz, Y = ln(-^Y (3) 

The appearance of this scale is a consequence of angular ordering (AO) of successive parton branchings 
in QCD cascades ISKl^. An important difference in the structure of light (i = q = u,d, s) versus heavy 
quark (h = Q = c,b) jets stems from the dynamical restriction on the phase space of primary gluon 
radiation in the heavy quark case, where the gluon radiation off an energetic quark Q with mass m and 



energy E ^ m is suppressed inside the forward cone with an opening angle 0^ = m/E, the above- 
mentioned dead cone phenomenon [4|. This effect is in close analogy to QED, where photon radiation is 
also suppressed at small angles with respect to a moving massive charged particle (e.g. tau versus muon). 

The corresponding evolution time parameter for a jet initiated by a heavy quark with energy E and mass 
771 appears in a natural way and reads [4] 

y = lii('^Y Kl = kl + z^m\ (4) 

which for collinear emissions © ^ 1 can also be rewritten in the form 



1 



K± = zQ, Q = E{e^ + el,)\ (5) 

with @ > @rn (see Fig.lTJ. 

An additional comment is in order concerning the AO for gluons emitted off the heavy quark. In ^, Q is 
the emission angle of the primary gluon g being emitted off the heavy quark. Now let Q' be the emission 
angle of a second gluon g' relative to the primary gluon with energy uj' <^ uj and ©" the emission angle 
relative to the heavy quark; in this case the incoherence condition ©'^ < (0^ + 0^) (see appendix lAb 
together with ©" > 0,„ (the emission angle of the second gluon should still be larger than the dead 
cone) naturally leads (|5]l to become the proper evolution parameter for the gluon subjet (for more details 
see H). For ©„, = 0, the standai^d AO (©' < 0) is recovered. Therefore, for a massless quark, the 
virtuality of the jet simply reduces to Q = EQ as given above. The same quantity k± determines the 
scale of the running coupling Ug in the gluon emission off the heavy quark. It can be related to the 
anomalous dimension of the process by 

7o^(K^) = 2iV,^iM = ^_^, /3o(n,) = -^f^iV,-|n,), A = ln-^, (6) 
IT /3o(y + A) 4iVc V 3 3 V Aqcd 

where nj is the number of active flavours and Nc the number of colours. The variation of the effective 
coupling Us as nj — )• nj + 1 over the heavy quarks threshold has been suggested by next-to-leading 
(NLO) calculations in the MS scheme [13] as well. However, since the inclusion of NLO terms to the 
couphng constant provides corrections of relative order 0{as) and therefore beyond C(^/a7), this effect 
is subleading in the MLLA (see ^). In this context /3o(?^/) will be evaluated at the total number of 
quarks we consider in our application. The four scales of the process are related as follows, 

(5 > m- > Qo ~ Aqcd, 

where Qq ~ ^qcd corresponds to the limiting spectrum approximation [51. Finally, the dead cone 
phenomenon imposes the following bounds of integration to the perturbative regime 

^<z<l-^, m'<Q'<E\el + el), (7) 

which now account for the phase-space of the heavy quark jet. The last inequality states that the minimal 
transverse momentum of the jet Q = EQm = m is given by the mass of the heavy quark, which enters 
the game as the natural cut-off parameter of the perturbative approach. 



3 Definitions and notation 

The multiplicity distribution is defined by the formula 

oo 
Pn = -F^oo = ) 2-,, -fn = 1 (8) 

where an denotes the cross section of an n-particle yield process, amei is the inelastic cross-section, and 
the sum runs over all possible values of n. 

It is often more convenient to represent multiplicity distributions by their moments. All such sets can be 
obtained from the generating functional Z{y, u) 111 defined by 

z{y,u) = Y,Pn{y){i + uT 

n=0 

at the energy scale y. For fixed y, we can drop this variable from the azimuthally averaged generating 
functional Z{u); the moments are then calculated from the MLLA master equation as 

1 d9^(e" - 1) 
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the average multiphcity is defined by the formula, 
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Fq are respectively the factorial moments, often called multiplicity correlators. Kg are the cumulants 
of rank q and Cq, the moments of the multiplicity distribution P„. The multiplicity correlator and the 
moment of rank g = 2 are related as follows, 

WA^-l))^ _ 

Moreover, the width D^ = (^N"^ ) — N"^ of the multiplicity distribution P„ can be written in the equivalent 
forms, 

In terms of Feynman diagrams, Fq correspond to the set of all graphs while Kq describe the connected 
diagrams. Therefore, Kq are more suited for the construction of the evolution equations. 

Specifically, we will compute the average multiplicity of partons in jets to be denoted hereafter as A^^. 
with A = Q,q, g, corresponding to a heavy, light quark or gluon initiated jet respectively. Likewise, we 
will compute the second rank multiplicity correlator inside the same jet. 

Once arrived at his point, let us make an important distinction between two different particle sources 
populating heavy-quark initiated jets. On the one hand, parton cascade from gluon emission yields the 
QCD component of the total jet multiplicity (the main object of our present study), excluding weak 
decay products of the leading primary quark at the final stage of hadronization. On the other hand, the 
latter products coming from the leading flavoured hadron should be taken into account in the measured 
multiplicities of jets. We shall denote the average charged hadron multiplicity from the latter source as 
A^^'^. Hence the total charged average multiplicity, A^j^*"', reads 

j^total ^ ^ch ^^f. A = q, Q. (9) 



As a consequence of the LPHD, N'^ = K,^^ x A^^ II12II13L where the free parameter IC^^ normahzes 
the average muhiplicity of partons to the average multiplicity of charged hadrons. For charm and bottom 
quarks, we will respectively set the values N^" = 2.60 ± 0.15 and N^"" = 5.55 ± 0.09 SIH]], while in 
hght quark jets one expects N^"" = 1.2 ± 0.1 ||T5l . 

Now let us point out the distinct trends from each contribution to © as the quark mass increases. The 
dead cone effect suppresses Nq for heavier quark masses. Conversely, N^ becomes more significant 
for bottom jets. As we shall later see, the former will ultimately dominate the behaviour of the total 
average multiplicity Nq*'^^ of heavy quark jets for high Q values. In this paper, we advocate the use of 
such a difference between average jet multiplicities as a signature to distinguish a posteriori heavy from 
light quark jets, particularly in fe-tagging techniques applied to the analysis of many interesting decay 
channels. 

4 QCD evolution equations 



The splitting functions fl^ 

P(z, as) = a,p(o) {z) + alP^^^ {z) + . . . (10) 

where P^^' (z) and P^^> {z) are respectively the LO and NLO splitting functions, can be associated to 
each vertex of the process in the partonic shower. P{z, as) determines the decay probability of a parent 
parton (quark, anti-quark, gluon) into two offspring partons of energy fractions z and 1 — z. In this 
paper, we are rather concerned with calculations which only involve the LO P^^^ {z) splitting functions 
in the evolution equations Q. At LO or tree level, the Lorentz structure of the massless and massive 
splitting functions is universal for the decays e^ — ;■ e^7, 7 — ;■ e+e~ and q ^ qg, g ^- qqin QED and 
QCD respectively ||9j[I21. The inclusion of the NLO splitting functions □ in our approach would provide 
corrections of relative order 0{as), which going beyond 0{^/a^) in ([T]), are subleading in this frame. 
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Figure 1: Parton splitting in the process Q — )• Qg: a dead cone with opening angle 0^ is schematically 
shown. 



Let us start by considering the the splitting process, Q — )• Qg, Q being a heavy quark and g the emitted 
gluon which is displayed in FigUl the corresponding LO splitting function reads lOlfTOl 
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Cf 
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z z{\ — z)m? 
2 fc^ + z'^m? 



P^Qli^) = P'oii^ 



(0), 



(11) 



''The dependence of our approach on a certain renormahzation scheme would arise with the inclusion of the NLO QCD 
splitting functions in the evolution. 



where k± ^ min{zEQ, (1 — z)EQ) is the transverse momentum of the soft gluon being emitted off the 
heavy quark. The previous formula (ITTI ) has the following physical interpretation, for k±_ <^ z^rr?, the 
corresponding limit reads PQg{z) — >• 2N'^ ^^'^ ^^^^ ^^ why, at leading logarithmic approximation, the 
forward emission of soft and coUinear gluons off the heavy quark becomes suppressed once Q <C &m^ 
while the emission of hard and coUinear gluons dominates in this region. 

For the massless process g — )■ gg, we adopt the standard three gluon vertex kernel 
and finally for g -^ QQ, we take llOlfTOl 



{l-z)[2-zil-z)], 



(12) 
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2z(l — z)m 
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(13) 



k'j_ + m^ 

which needs to be resummed together with the three gluon vertex contribution. However, as a first 
approach to this problem, we neglect the production of heavy quark pairs inside gluon and quark jets, 
making use of | 



p^'Az) = p'M 
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rra=0 



[1 - 2z(l - z)] 



(14) 



gq V-/ ~ ^ gQ K-ni'i—u ^^ 

Including mass effects in the evolution equations also requires the replacement of the massless quark 
propagator l/k'j_ by the massive quark propagator l/{k'^ + z'^m?) fTO'l, such that the phase space for 
soft and coUinear gluon emissions off the heavy quark can be written in the form lITOl 



dVgg 



70^ 



dkl 



k^ + 
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^dzP^Qli-l 
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where PX (z) is given by (fTTI) . Working out the structure of (ITTI) and setting k^ ss zEQ one has 
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such that, one can recover the phase space for the soft and coUinear gluon emission in the DLA 
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Notice that the overall multiplicity of the process e'^e'^ — ?■ QQg cannot be represented simply as the 
sum of three independent parton multiplicities lU. As stressed in Q, the accompanying multiplicity off 
the quark dipole becomes dependent on the geometry of the whole jet ensemble in a Lorentz invariant 
way and should be treated as a different problem. However, as suggested in Q and demonstrated in the 
appendix |BJ we recover the correct limit of the one jet event ( 115116! ) from the heavy quark dipole |4|, in 
the soft z <^ 1 and coUinear 9^1 gluon limits. 

According to the Low-Bumett-Kroll theorem lil9[|20L the dz/z part of the radiation density has a clas- 
sical origin and is, therefore, universal, independent of the intrinsic quantum numbers and the process, 
while the other terms are quantum corrections. The system of two-coupled evolution equations for the 
gluon and quark jets average multiplicity in the massless case at MLLA simplifies to the following in the 
hard splitting region fc^ ~ £'9 (z ~ 1 - z ~ 1) 1118112111 
d fl-Qo/Q 

dY' 



:NqiY) 



_d_ 
'dY 



Ng{Y) 



Qo/Q 

i-Qo/Q 
Qo/Q 



dzji{z)p(^\z)Ng{Y + lnz) 



(18) 



dzjiiz) fpW(z)iV,(y + lnz)+ nM^Xz) {2N,{Y) - Ng{Y)) 



,(19) 



where Pgg' {z) = Pq (2;)|m=o and Pgg' {z) = PJ^J {z)\m=o according to (fTTI ) and ( |T3l ) respectively. It 
is obtained from the MLLA master equation for the azimuthally averaged generating functional Z{y, u), 
by taking the functional derivative over u (see section |3]l. The arguments of A^^ in the right hand side 
of the equations do not depend on z because for hard partons z ~ 1, the original arguments y + Inz 
and Y + ln(l — z) of these functions can be approximated to Y after In z and ln(l — z) are neglected. 
Substituting (fT6l ) into (fTSl ). replacing EQ by E{Q'^ + 6^)^^^ in the argument of all logs, which is 
equivalent to replacing the massless propagator by the massive one and the argument of the running 
coupling k±_ = zEQ by k± = zE{Q^ + 0^)^'^, after taking the bounds (|7]) and integrating over the 
regular part of the splitting function ( fT6l ). one has the equation 
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finally, introducing the chain of transformations 
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-2Y+2Lr, 



Lir, 



, m 



where the Lm logarithms are new in this context and provide power suppressed corrections to the solution 
of the evolution equations. Since massless, in the gluon jet however, the evolution time variable remains 
the same Q as that in the massless quark jet. Nevertheless, the argument of the average multiplicity in 
the gluon subjet Ng{y) in (l20l ). depends on the same argument y. That is why, in the following, we insert 
the mass of the heavy quark in the argument of all logs in (fT9l) and take the same integration bounds ([T]) 
correspondingly. The system of QCD evolution equations now reads 
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(21) 
(22) 

(23) 

(24) 
(25) 



There are the following kinds of power suppressed corrections to the heavy quark multiplicity: the lead- 
ing integral term of (|2TI ) is 0{^) suppressed, while subleading MLLA corrections appear in the stan- 



dard form 0{y/a2) like in the massless case, finally O(^y^) and O(^y^), which are new in this 
context. Similar corrections have been found in the treatment of multiparticle production off the heavy 
quark dipole 15] and in the computation of the heavy quark content inside gluon jets [22]. 



Similar power suppressed corrections proportional to 6^^+^™ in A{Y,Lm) and e2{Y,Lm) were re- 
ported in [21 1 for the massless case. Indeed, such results can be recovered after setting m/Q — )• Qo/Q 
in ( |23] ) and (|25] ). For massive particles however, these terms are somewhat larger and can not be ne- 
glected in our approach unless they are evaluated for much higher energies than at present colliders. On 
top of that, the corresponding massless equations in the high energy limit are obtained from (|2TI) and 

~ ~ ^ y, y„ ^ 0, ei ^ 1, 



(1221) simply by 


setting y-^y,Y -^Y, Y^v 
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62 ^ 62 = - - 


. 3^-^ + O (e-2^) , A^ 


and are written in the standard form 1 18112 11 



A = airii 



2 + ^fl-2^^ 



2N,, \ Nr, 



g-y , rn (.~2Y-, 



+ O (e-2^) , 



^^ = 1 dy^l{y)Ng{y)-UYhl{Y)N,{Y), (26) 

^ = j\y^l{y)N,{y)-A{Y)^l{Y)N,{Yl (27) 

with the initial condition Ng^q{Y = 0) = 1 at threshold. Notice that ((2TI) and ((22l) are valid only for 
m^ Qq and therefore ?7i — )• does not reproduce the correct limit, which has to be smooth as given by 
the massless equations (l26l) and (l27l) . 

As can be seen from (|2TI) . the function ei also gives the power suppressed contribution oc _e~2y+2Lm 
which decreases the production of soft and coUinear gluons off the heavy quark, however, this contri- 

2 

bution is power suppressed 0{^) and turns out to be rather small as the energy scale increases. Since 
heavy quarks are less sensitive to recoil effects, the subtraction terms oc e"^"*"^™ and oc e-'2Y+2Lm ^^ 
e2(Y, Lm) diminish the role of energy conservation as compared to massless quark jets. As a consistency 
check, upon integration over Y of the DLA term in Eq.(l2TI). the phase space structure of the radiated 
quanta in dTT] ) is recovered: 

:i2 






Im/Q 

Notice that the lower bound over 6^ in (|28]) (Y in ^T^) can be taken down to "0" {Ym = Lm in (|2T]) ) 
because the heavy quark mass plays the role of collinear cut-off parameter. 



4.1 Towards the solution of the evolution equations 



First we solve the self-contained equation for the gluon jet (122] ). The second and third exponential terms 
in (|23] ) are slowly varying functions of the variable Y at high energy scales. In the same Umit and for 
the sake of simplicity we set the bound of integration over {z,y) to Y^v — ?• Y , Y^ — )• and solve the 
equation by performing the Mellin transform 

Ng(Y) = I ^e^yNgioj), (29) 

where the contour C lies to the right of all singularities in the complex plane of uj. Replacing (l29l ) into 
((22]) leads the first order differential equation in Mellin space: 



which upon inversion leads to the following solution 



Ng{Y,Lm) = const x{Y + Xy^exp 2W^^ . (31) 

Then, the initial condition at threshold, which is reached when the jet virtuality approaches the mass of 
the heavy quark yields 



\S_„ I o, /^m + '^ 



Ng{L^) = 1 =^ const = (y„ + A)^ exp | -2W ^ | (32) 



and finally, 



-s 



with 



From (l33l) and by making use of Ng (Y) ~ exp (J dy^{y)\ |[5l. one gets the rate of multiplicity growth 
as a function of Y to be, 

7 ~ 70 - ^ll- 

A similar solution without power corrections, which was written for the "gluon mass" was given in |[22l . 
Notice that the fact of introducing the "gluon mass" in this context is technical rather than physical. 
However, phenomenological observations favour a dynamically generated mass for the gluon |[23]| . In 
order to obtain the approximate solution of (|2TI ). as before we consider the functions 

as constants at high energy scale. Subtracting (|2TI ) from (l22l ) and setting 

N, d ( idNQ\ o ^ 

on the right hand side of the outcoming, one has 

dY Cf ^ dY Cf^ ^ ^ dY \^ dY J 

Working out the structure of (l34l ). so as to obtain the ratio 



we can rewrite it in the form 
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Finally, we obtain 

r{Y,Lm) ^^' ^e^\Y,L„,) [l - (A(y,L„) - ei\Y,L^)e2{Y,Lrn)) 70 



(36) 



(37) 



which becomes vahd in the limit 

^ ^^' 0, n{Y,Lm) "^^^ (a{Y,L^) - ei\Y,K,)e2{Y,Lm)) . (38) 

dY ^ ' 

Finally, the approximate average multiplicity in a jet initiated by a heavy quark reads 

iVQ(^,^™) = ^^i(^'^™)7^^^^^^^' (39) 

^^c 1 -ri(y,L„,)7o 

where ei(y, L^) is written in ( [251 ) and Ng(Y, L^) is given by ( [33l ). Thus, as the mass of the leading 
heavy quark increases, the multiparticle yield in the heavy quark jet is affected by power corrections, 
by the suppression of the anomalous dimension 70 = 70 (m^) and mainly by the massive suppressed 
exponential contribution arising from the initial condition at threshold. However, for the sake of com- 
pleteness, we solve the evolution equations numerically and display the energy dependence of the average 
multiplicity in Figj2] The asymptotic behaviour of the distribution is then seen to follow the expected 
exponential increase given by ( [39l ). with Ng in (|33] ). Finally, we estimate the difference between the light 
and heavy quark jet multiplicities, which yields. 



Ng-NQ « 




1 — exp — 2W — — 



Ng, Ng(xexp2\ —. (40) 



Hence, (1401 ) is exponentially increasing because it is dominated by the leading DLA energy dependence 
of Ng. According to (l40l ). the gap arising from the dead cone effect should be bigger for the b than for the 
c quark at the primary state bremsstrahlung radiation off the heavy quark jet. The approximated solution 
of the evolution equations leads to the rough behaviour of Ng — Nq in (l40l ). which is not exact in its 
present form. In FiglJl we display the numerical solution of the evolution equations (l26l ) for Ng and (12T1) 
for Nq corresponding to the heavy quark mass intervals ttt-c = 1 — 1.5 GeV, m^ = 3 — 5 GeV. Let us 
remark that the gap arising between the light quark jet multiplicity and the heavy quark jet multiplicity 
follows the trends given by (l40l ) asymptotically with £' — ;■ 00. In particular, notice that the dispersion 
of the mean multiplicities becomes irrelevant for the purposes of our study. This behaviour should not 
be confused with that followed by the QQ antenna in the e+e^ annihilation, where the difference is 
roughly constant and energy independent ||4l|24l. Indeed, (l40l ) can not be extrapolated to the dipole case 
by simply setting Nq q w 2Nq because the evolution equations do not take into account interference 
effects between the Q and the Q jets in the e+e^ annihilation. Finally, as expected for massless quarks 
Lm = 0, the difference Ng — Nq vanishes. 

5 Heavy quark evolution of second multiplicity correlator 

The second multiplicity correlator was first considered for massless quarks in [25]. It is defined in the 

(2) 
form N^ = {Na{Na — 1)) in gluon (A = g) and quark {A = q) jets. The normalized second 

multiplicity correlator defines the width of the multiplicity distribution and is related to its dispersion 

squared D\ = (Nj^) — N'^ by the formula (see definitions and notation in section[3l) 

D\ = {Fa,2 - ml + Na. (41) 
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Figure 2: Massless and massive quark jet average multiplicity Nq as a function of the jet hardness Q. 
Bands indicate nic and ?ti{, in the [1, 1.5] and [3, 5] GeV intervals respectively. 



The second multiplicity correlators normalized to their own average multiplicity squared are 



^2,9 



G, 



{NgiN, - 1)) 



^1 



'2,q 



Qi 



{N,iN, - 1)) 



Ar2 

9 



(42) 



inside a gluon and a quark jet respectively. These observables are obtained by integrating the double 
differential inclusive cross section over the energy fractions xi = ei/E and X2 = 62/ E of two particles 
emitted inside the jet, 

{Na{Na - 1)) = [[dxidx2 (-- '^^'^ 



^a dxidx2 J A 

The system of evolution equations for light quarks following from the MLLA master equation can be 
written as ll5l[T8l. 



A (atW _ Ar2) = r ^"/^ dz^l{z)PJS^ iz)Np (Y + In .), 
"^ JQo/Q 

A (Ar(2) _ Ar2) = f'"'''^'' dzjliz)p('^^ {z)Nf) {Y + In .) 
"^ JQq/Q 



(43) 



+ nf 



l-Qo/Q 



Oo/Q 



dz7o' (^)^i? (^) [2 (<) (y ) - Ar2 (y )) _ (atC^) (y ) - iV^ (y ) 



+ (2iv,(y)-iv,(y) 



(44) 



with the following relations at DLA ||26[|27i . 



Ar(2) _ Ar2 = ^ fAr(2) _ n^\ m =—N„ 



r(2) 



(45) 



The arguments of N)^ and Nj^ in the right hand side of the equations do not depend on z because 
for hard partons 2; ~ 1, the original arguments y + Inz and Y + ln(l — z) of these functions can be 
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approximated to Y after In z and ln(l — z) are neglected. Substituting ( fT6l ) into (|43] ). after replacing £^0 
by E{Q'^ + 0^)^/^ in the argument of all logs, taking the bounds Q and integrating over the regular 
part of the splitting functions, one has the system 



dY 



(iV(2) _ AT^ 



y^v 



irn 



dy^l{Y)N^'){y) - A{Y ,L„H(X)N^^\Y) 



+ [A{Y,L^) - B{Y,L,^)hl{Y)Nl{Y), 



(47) 



with the initial conditions N\ ' (L. 
B{Y,Lm) = binf) 



dN' 



(2) 



a^iLni) =0, where 



2 J- 1-2 ^ 



with 



2Nr 



b{nf) 



Nr 



^-Y+L^ + 



l-^(l-2^^ 



A^. 



Nr 



-2Y+2Lr, 



AN, 



11 2 ( Cp 



Accordingly, in the massless limit, (I46b and (1471) reduce to 



(48) 



CpdY' 



Y 



(iV(2) - Ar2) = / dy^l{y)N^^\y) - UY)iI{Y)N^'\Y), 



A(Ar(2) _ Nl) = £ dy^l{y)Nf\y) - A{Y)^l{Y)Nf\Y) 



+ [A{Y)-B{Y)hl{Y)Nl{Y), 



(49) 



(50) 



where 



B{Y) = b{nf) 



nr I Cp 

2 ^ 1-2 



2Nr 



Nr 



+ 0(e-^^). 



The functions £2(1^) and AiY) are defined above through the equations for the average multiplicity ([26 
and ( |27l ) in light quark jets. 



5.1 Approximate solution of the evolution equations 



For the gluon jet, taking into account that at high energy scales one has A w A{Y ^ Lm) and B ?» 
B{Y, Lm) (^^ ^^^ 0), and making use of ^, the solution reads lH 



G2{Y,Lm) - 1 ""^^ ^ - Ci(y,L„)7o, 



(51) 



where. 



Accordingly, for the quark jet one finds 129 



Ci(y,L„) ^^^ -|A(y,L„) + i/3o(n;) + \B{Y,Lm), ^ ^^^ 0. 



Q2{Y,Lr^ 



,%.§.,.,>■,.„.) Q- 



C'i(y,L^)7o 



(52) 
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where 






+ ^A(n/) + iB(F.Z,„; ''^'"-=*' 



0. 



are mainly affected by power corrections 0{^y/a^) and 



Therefore, the correlators G2 ( |5T| ) and (52 

2 
0(^1/0!) which diminish the role of energy conservation in a heavy quark jet and make the correlation 

stronger as the particle yield gets suppressed inside the dead cone region. Thus, with such effects, the 
correlators increase as the mass of the leading heavy quark increases and approach the asymptotic DLA 
values G2 = I and Q2 = ^ + 3^ respectively. However, for realistic energy scales this approximation 
fails, in particular because of the integration over the dead cone term oc ei{Y, Lm) in the double loga- 
rithmic contribution of (l46l ). That is why, as we further emphasize in the appendix O we should rather 
display the numerical solution of the equations (|46l ) and (|47] ) in the relevant energy range. 



6 Phenomenological consequences 

The study of multiplicity distributions (mean and higher rank momenta) and inclusive correlations has 
been traditionally employed in the analysis of multiparticle production in high energy hadron collisions, 
notably regarding soft (low pt) physics (see e.g. ifTSi and references therein). Moreover, the use of 
inclusive particle correlations has been recently advocated in the search of new phenomena ||30l . 

On the other hand, it is well known that the study of average charged hadron multiplicities of jets in 
e+e" collisions has also become a useful tool for testing (perturbative) QCD calculations (see B1I241 
and references therein). 




Figure 3: Massless and massive quark jet average multiplicity _/Vq*"' as a function of the jet hardness Q 
including heavy quark flavour decays. Same comments as in Fig.|2l 



In this paper we advocate the role of mean multiplicities of jets as a potentially useful signature of new 
physics when combined with other selection criteria. In Fig.O we plot as function of the jet hardness Qj, 
the total average jet multiplicity Q, which accounts for the primary state radiation off the heavy quark 



'The energy range 100 < Q{GeY) < 200 should be realistic for Tevatron and LHC phenomenology. 
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together with the decay products from the final-state flavoured hadrons, which were introduced in section 
[3] For these predictions, we set K.^^ = 0.6 in ©, which we take from 1,31,1 and Qo ~ ^QCD = 230 MeV 
ll5l . Moreover, the flavour decays constants N^'^ = 2.60±0.15 and Nj^"^ = 5.55ib0.09 are independent of 
the hard process inside the cascade, such that N'^ can be added in the whole energy range. For instance, 
such values were obtained by the OPAL collaboration at the Z^ peak of the e+e^ annihilation. In this 
experiment, D* mesons were properly reconstructed in order to provide samples of events with varying c 
and b purity. By studying the charged hadron multiplicity in conjunction with samples of varying b purity, 
it became possible to measure light and heavy quark charged hadron multiplicities separately |[T4l . As 
compared to the average multiplicities of the primary state radiation displayed in Fig.[2j after accounting 
for N'^, the b quark jet multiplicity becomes slightly higher than the c quark jet multiplicity, although 
both remain suppressed because of the dead cone effect. 




20 40 60 



100 120 140 160 180 200 

Q (GeV) 



Figure 4: Massless and massive quark jet correlator Q2 as a function of the jet hardness Q. Same 
comments as in Fig.|2l 



The second quark jet correlator defined in (l42l ) for different flavours is displayed in Fig.|4]as a function 
of the jet hardness Q. Since contributions to the dispersion from quark flavour decays are negligible 
(A'^'^ = 2.60 lb 0.15 and A'^^'^ = 5.55 ± 0.09) the correlation is the strongest for partons at the primary 
state radiation of the process. Notice that, while the u, d, s and the c quark correlators are of the same 
order of magnitude for a jet hardness Q ^ 40 GeV as relevant energy range, the vertical difference 
with the b quark correlator, which is weaker, can still exceed ~ 20%. Therefore, the measurement of the 
quark correlator should provide a further signature of b flavour and associated exotic particles yield when 
compared with u, d, s, c correlators. Finally, the variation of the charged hadrons average multiplicities 
and the correlator in the above-mentioned intervals for the charm and and bottom masses turns out to be 
neghgible for our purposes. 
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7 Conclusions 

Jet physics has been so far of paramount importance in the rise and development of the SM and ex- 
pectedly will keep such a prominent role in the discovery of new phenomena at hadron colliders like 
the Tevatron and the LHC. However, QCD jets represent a formidable challenge to disentangle signals of 
new physics from hadronic background in most cases. On the other hand, plenty of new physics channels 
end with heavy flavours in the final state, before fragmenting and hadronizing. 

Thus, our present work focusing on the differences of the average charged hadron multiplicity between 
jets initiated by gluons, light or heavy quarks could indeed represent a helpful auxiliary criterion to 
tag such heavy flavours from background for jet hardness Q > 40 GeV. Notice that we are suggesting 
as a potential signature the a posteriori comparison between average jet multiplicities corresponding 
to different samples of events where other criteria to discriminate heavy from light quark initiated jets 
were first applied. In other words, one should compare mean multiplicities at different jet-hardness Q, 
in order to check that they agree with QCD predictions. FigjS] plainly demonstrate that the separation 
between light quark jets and heavy quark jets is allowed above a few tens of GeV with the foreseen 
errors of the experimentally measured average multiplicities of jets. The difference between light quark 
jet multiplicities and heavy quark jet multiplicities Nq—Nq in one jet is exponentially increasing because 
of suppression of forward gluons in the angular region around the heavy quark direction. This result is 
not drastically affected after accounting for heavy flavour decays multiplicities, such that it can still 
be used as an important signature for the search of new physics in a jet together with other selection 
criteria. However, our result can only be applied to single jets and therefore, it should not be extrapolated 
to the phenomenology of the QQ dipole treated in H because neither interference effects with other 
jets nor large angle gluon emissions are considered in our case. As a complementary observable, in 
particular for 6-tagging, the second multiplicity correlator (l42l) displayed in Fig.|4] should also contribute 
to discriminate h quark from u, d, s, c quark channels. Indeed, while the c quark correlator remains 
of the same order of magnitude than the light quark jet correlator, the b quark correlator gets weaker 
by 20% and therefore, distinguishable with respect to the other quarks in the relevant energy range. 
Furthermore, the inclusion of the heavy quark mass in the evolution equations for the correlator does not 
affect the asymptotic energy independent flattening of the slope arising from the KNO (Koba-Nielsen- 
Olsen) scaling B2I . 

Notice that the measurement of such observables require the previous reconstruction of jets at hadron 
colliders. Thanks to important recent developments on jet reconstruction algorithms |[33l435l . future 
analysis such as single inclusive hadron production inside light and heavy quark jets look very promising. 
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A From AO to the incoherence condition of gluon emission off the heavy 
quark 



In the MLLA, the parton decay probabihties are written in a form m, 



dn 



dl^^<faQ^Q, = ^P^l{z)dzV{n)'^. VI, 



Svr- 



?(Q) 



.n 



flg'Q + agQ - O-g'g! 



(53) 





^' 


/ ^'^^^ 


/ 
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1 


Q=(c, b) "" - - ~ ^ ^ ^ ^ \ 



Figure 5: Second gluon emission off the primary gluon in the process Q — ;• Q + 5(/) + g\s). 



It describes the process Q — )■ Q + 5(/) + ^'(s) displayed in FigJH where the subscripts mean father 
and son. In this case we define = O^q, 0' = ©g'g- For light quarks involved in the same process 
g — )■ g + g'(/) + g\s), if "i" and "k" denote the massless particles, then the angular factor a^^ in the 
relativistic case is written as 

fljfc = 1 -cos6ifc. (54) 

After taking the azimuthal average around the "son" gluon direction, one obtains lU, 

.2^ AA. . 2 



dig) 






27r ^'■•^i aai 



■diflgq -agig), 



(55) 



a a 



with "9 the Heaviside function. This leads to the exact AO inside partonic cascades by replacing the 
strong AO in the DLA 0' <C by 0' < in the MLLA. For massive particles we may write (1531 ) in the 
same form after replacing the standard massless splitting functions [SJ by the massive one ifTOll . If the 
leading parton is a heavy quark, the angular factor of the emitted gluon "(7" off the heavy quark Q, can 
be checked after some simple kinematics, to be written in the form. 



OgQ = 1 - \/l - 0-m cos 0, 
where 0m is the angle of the dead cone. In this case, ( |55l) can be rewritten as follows 

2 



V^ 

9(Q) 



27r 



^''aiQ)^ 



a„i 



--^iagQ - ttg'g), 



(56) 



(57) 



9 9 



imposing that cos 0' > y^l — 0^ cos 0. For small angles, if one sets cos Ri 1 — ^ in both members 
of the previous inequality, one gets the incoherent condition: 



)2. 



e'2 < 02 + e^ 

In the massless case 0^ = 0, ( [58l ) simply reduces to the standard exact AO 0' < 0. 



(58) 
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B Accompanying radiated quanta off the heavy quark dipole 



In milBl, the probability of soft gluon emission of the heavy quark pair QQ produced in the e~^e 
annihilation was written in the form, 



d^a, 



QQg 



CpOisiK^) dz /3 



TT 



Z V 



—dcos ©c < 2(1 — z 



13^ sin2 e. 



(l-/32cOs2ee)2 



+ z' 



1 



1 



1 - /32 C0S2 e^ 



Cv ( ' 
(59) 

where z is the energy fraction of the emitted gluon and Gc the emission angle with respect to the center 
of mass of the QQ pair. Moreover, the following notation was introduced: 

4m2 



/3^ 



1 



v^ = l 



4m2 m? 



(60) 



W^il-z)' 

The transverse momentum of the gluon appearing on the argument of the running coupling in (|59l l was 
written in the form. 



(61) 



^) (l-/3^cos2G,)^ 

With such a notation, we now take interest in the soft (1 — z ~ 1) and coUinear (Be ^ 1) limits of ( [59l ) 
and set W — )• 2E in order to reduce ( [59l ) to the single jet event initiated by a heavy quark Q. Thus, the 
terms take the following form 



P^'^' 



m 



1 - el, (1 - p' cos^ G,)2 ^<<'^^<<' (el + elf, 



2(1 - z)/32 sin2 Gc 2«i^c«i 2(1 - z)e2 

(i-/32cos2ee)2 ^ Jef+eif' 



1 - /32 C0S2 Oc 



a 



2<i,ec<i 



92 + 92 ■ 



The term proportional to — ^ in the cross section (l59l ) does contribute neither as a soft logarithmic contri- 
bution nor as a collinear one, and therefore can be neglected in this approximation. It should correspond 
to a Feynman diagram which only accounts for interference effects between the Q and the Q jets in the 
QQ antenna. In this limit, for one jet we set d'^aqq -^ d'^aqg, 9c — ;> 9 and ^ — > 1, such that the 
dipole cross section ( [59l ) can be rewritten in the form 

d aqg ~ irr7o(«J 



Nr. 



92 



Z 

92 + 92 1 z 92 + 92 " + 2 ^ 92 + 92 



dz 



(62) 



- ll{n'^^-^P^l{z)dz, 
1% 



M 



as given in (IT5] ). where Pi (z) was written in (IT6l ). while (|6TI ) was reduced to the following, 



Therefore, in the soft and collinear approximation, the dipole case ( 1591 ) reduces to the jet event (1621) . 
which coincides with the expression given in (ITSl) . In the massless case 9m = 0, it simplifies to the 
standard DGLAP kernel as explained in |4|, 
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C Analytical versus numerical solution of the heavy quark correlator 
equation (1461) 

In Fig.|6]we display the analytical solution together with the numerical solution of (l46l ) for the second 
multiplicity correlator Q2 defined in (l42l ). As it can be seen, when the mass of the leading heavy quark 
increases, the approximated analytical correlator becomes slightly stronger. However, because of forward 
gluon suppression taken into account by the integrated function ei {Y , Lm) in the leading DL contribution 
of (l46l ). such a behaviour cannot be trusted for lower virtualities than few thousands of GeV. That is why, 
even if the shape of the analytical solution may be correct for Q > 100 GeV, we should only trust the 
shape and normalization of the numerical solution in a much wider energy range Q > 40 GeV in view 
of realistic QCD predictions. 
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Figure 6: Analytical (left) versus numerical (right) solution of equation (1461 ) for the second multiplicity 
correlator Q2 defined in (l42b . 
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